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Abstract 

The two neutrino double beta decay of 124 ' 126 Xe, 128 ' 130 Te, 130 ' 132 Ba and 150 Nd isotopes is studied in 
the Projected Hartree-Fock-Bogoliubov (PHFB) model. Theoretical 2v (3~ (3~ half-lives of 128 . 130 Te, and 
150 Nd isotopes, and 2v (3+ (3+ , 2v (3+ EC and 2v EC EC for 124 ' 126 Xe and 130 ' 132 Ba nuclei are presented. 
Calculated quadrupolar transition probabilities B(E2: + — > 2 + ), static quadrupole moments and g factors 
in the parent and daughter nuclei reproduce the experimental information, validating the reliability of the 
model wave functions. The anticorrelation between nuclear deformation and the nuclear transition matrix 
element is confirmed. 

PACS Numbers: 23.40.Hc, 21.60.Jz, 23.20.-g, 27.60.+j 

Introduction 



nuclear double beta decay ((3(3) is a weak process in which two neutrons (protons) inside a nucleus decay 
j^into two protons (neutrons) , emitting two electrons (positrons) [TJ [2] . In the lepton- number-conserving process, 
f^*namely two neutrino double beta (2v (3(3) decay two antineutrinos (neutrinos) are also emitted [3-8]. It has 
^been observed in several nuclei [10-13]. Given that lepton number is not a gauge symmetry, several extensions 
. r^oi the standard model predict its violation. In this neutrinoless double beta (Of (30) decay no neutrinos are 
^^rnitted. This Qv (3(3 decay has never been observed, although controversial evidence has been published [14] . 
'"jl positively detected, it would provide unique information about the Majorana mass of the neutrino, and for 
■ - the absolute scale of neutrino masses [15-20]. 

Theoretical calculations of (3(3 decay represent one of the hardest challenges ever faced in nuclear physics. 
The process is definitively non-collective, and proceed through strongly suppressed channels, which are very 
sensitive to details of the wave functions of the parent and daughter nuclei. While specific models have been built 
which are able to adjust (postdict) most of the observed 2v (3(3 decay half lives, they usually have parameters 
which allow different results as well, making their predictions unreliable. A case by case analysis, studying each 
nuclei with the best available model, and describing as much nuclear observables as possible within the same 
theoretical scheme, seems to be a sensible approach [21-23]. 

The 2v (3(3 decay can occur in four different modes, namely double electron ((3~ (3~) emission, double positron 
((3 + (3 + ) emission, electron-positron conversion ((3 + EC) and double electron capture (EC EC). The later three 
processes are energetically competing and we shall refer to them all generically as positron double beta decay 
(e + DBD) modes. Nuclear matrix elements (NMEs) Mi v associated with the 2v (3(3 decay can be extracted 
directly from the observed half-lives of eleven nuclei undergoing 2v(3~ (3~ decay out of 35 possible candidates 
[T3"] . In case of 2v e + DBD modes, experimental limits on half-lives have already been given for 24 out of 34 
possible isotopes [13]. 

In all (3(3 decay emitters, which are even-Z and even- TV nuclei, the pairing of like nucleons plays a funda- 
mental role, energetically inhibiting the decay to the intermediate odd-odd nuclei. The quadrupole-quadrupole 
interaction drives nuclei to deformed shapes when both proton and neutron shells are open. The interplay be- 
tween this two major components of the effective nuclear interaction can generate quite complex energy spectra. 
In the mass region Aml30 it can be seen that, while Te isotopes have a vibrational excitation spectra, Xe and 



Ba isotopes develop rotational bands. The mass region A ~ 150 offers an example of shape transition i.e. the 
sudden onset of deformation at neutron number N—90. Nuclei range from spherical to well deformed, with large 
static quadrupole moments. Thus, it is expected that pairing and deformation degrees of freedom will play some 
crucial role in the structure of i24,i26,i28,i30,i32 Xej 124,120,128,130^ 130,132^ i 5 o Nd and i50 Sm nuclei xhe p 

decay can be studied in the same framework as many other nuclear properties and decays. A vast amount of 
data has been collected over the past years concerning the level energies as well as electromagnetic properties 
through experimental studies involving in-beam 7-ray spectroscopy. The availability of data permits a rigorous 
and detailed critique of the ingredients of the microscopic model that seeks to provide a description of nuclear 
(3(3 decay. 

Theoretical studies predict that deformation plays a crucial role in case of 2v (3~ (3~ decay of 100 Mo and 
150 Nd [231 US]. Auerbach et al. [SB] and Troltenier et al. [27J have already shown that there exists an inverse 
correlation between the Gamow- Teller (GT) strength and quadrupole moment. The effect of deformation on 
the distribution of the Gamow- Teller and /3-decay properties has been studied using a quasiparticle Tamm- 
Dancoff approximation (TDA) based on deformed Hartree-Fock (DHF) calculations with Skyrme interactions 
[28] . a deformed self consistent HF+RPA method with Skyrme type interactions [29 . The comparison of 
the experimental GT strength distribution B(GT) from its decay with the results of QRPA calculations was 
employed as a novel method of deducing the deformation of the N = Z nucleus 76 Sr [3U]. The effect of 
deformation on the 2v (3~ f}~ decay for ground-state transition 76 Ge — > 76 Se was studied in the framework of 
the deformed QRPA with separable GT residual interaction [51] . A deformed QRPA formalism to describe 
simultaneously the energy distributions of the single beta GT strength and the 2v (3(3 decay matrix elements, 
using deformed Woods-Saxon potentials and deformed Skyrme Hartree-Fock mean fields was developed [32] ■ 

In all these works calculations are performed in the intrinsic basis, where angular momentum is not a good 
quantum number. The Projected Hartree-Fock-Bogoliubov (PHFB) model offers, in this sense, a competitive 
alternative. On one hand, the PHFB model restores the rotational symmetry, providing very reliable wave 
functions for the parent and daughter (3(3 decay emitters. On the other hand, in its present version the PHFB 
model is unable to provide information about the structure of the intermediate odd-odd nuclei, and, in particular, 
on the single (3 decay rates and the distribution of Gamow- Teller strength. Notwithstanding, the PHFB model 
has been successfully applied to the 2v (3~(3~ decay of many emitters in the mass region A ~ 100, where it was 
possible to describe, in the same context, the lowest excited states of the parent and daughter nuclei, as well as 
their electromagnetic transition strengths, and to reproduce their measured (3(3 decay rates on the other [55] , 

A quantitative description of the 2v (3~ (3~ decay [33J and 2v e + DBD modes [33] of nuclei in the mass region 
A ~100 for the + — » + transition has been obtained. The same formalism allowed the analysis of other 
observed nuclear properties, including the yrast spectra, the reduced transition probabilities B(E2: + — » 2 + ), 
the static quadrupole moments Q(2 + ) and the g -factors g(2 + ) of both parent and daughter nuclei. The 
study was performed using the PHFB model. Its application, in conjunction with the summation method, has 
motivated the present study of the 2v (3~ /3~ decay of 128 < 130 Te, 150 Nd and 2v c+DBD modes of 124 < 126 Xe and 
I30,i32g a isotopes. The 2v e + DBD modes of 124 ' 126 Xe and 130 ' 132 Ba isotopes were studied in the earlier work 
of Shukla et al. [35] . In the present work, the HFB wave functions are generated with higher accuracy and in 
these nuclei, it is observed that the results are very much dependent on the later. 

By the use of the pairing plus quadrupole-quadrupole PPQQ interaction ,36] the interplay between sphericity 
and deformation can be studied. In this way, the PHFB formalism, employed in conjunction with the PPQQ 
interaction, is a convenient choice to examine the explicit role of deformation on the NTME M 2v . A strong 
dependence of the M 2v on the quadrupole deformation was found varying the strength of the quadrupole- 
quadrupole interaction, for the 2v (3~ (3~ decay [33] and 2v e + DBD modes [34] of nuclei in the mass region 
A ^100. The anticorrelation between the quadrupolar deformation and the (3(3 decay transition amplitude is 
analyzed in detail. 

The present paper is organized as follows. The theoretical formalism to calculate the half-lives of 2v (3(3 
modes has been given in a number of reviews [121 HE| • Details of the PHFB study were presented in previous 
publications, for 2v (3~ ' (3~ decay [33] and 2v e + DBD modes [33] of nuclei in the mass region A ~100. Details 
of the mathematical expressions used to calculate the spectroscopic properties of nuclei in the PHFB model 
have been given by Dixit et al. [37]. Here, we briefly outline the relevant results in Section 2. In Section 3, 
we present a detailed study of the wave functions of i24,i26,i28,i30,i32 Xe) 124,120,128,130^ i30,i32 B&; i 5 o Nd and 

150 Sm nuclei, calculating the yrast spectra, reduced B(E2: + — > 2 + ) transition probabilities, static quadrupole 
moments Q{2 + ) and g -factors g{2 + ) and comparing them with the available experimental data. The half-lives 



T*J 2 for the 2v (3" (3" decay of 128 < 130 Te, 150 Nd and 2v e+ DBD modes of 124 ' 126 Xe and 130 < 132 Ba isotopes for 
the + -^0 + transition are calculated. The role of deformation on NTMEs M 2v is also studied. We present 
some concluding remarks in Section 4. 

2 Theoretical framework 

The inverse half-life of the 2v (3(3 decay modes for the + — > + transition is given by 



= G 2v {k)\M 2v (k)f 



(1) 



where k denotes the (3 (3 and e + DBD modes. The integrated kinematical factor G 2v (k) can be calculated 
with good accuracy [T5] and the NTME M 2v (k) is defined as 
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where the total released energy Wo(k) — 2Eq, for the different 2v (3(3 decay modes, is given by 
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(3) 
(4) 
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The numerical evaluation of Eq. ((2|) requires an explicit summation over all the states of the intermediate odd 
Z-odd N nuclei. However, it is not possible to study the structure of intermediate odd-odd nuclei in the present 
version of the PHFB model. Therefore, we carry out the summation over intermediate states [33l [34] using 
the summation method given by Civitarese et al. [35] • This technique is similar to the Operator Expansion 
Method [35], which has been proven to be in general inaccurate [35J SDJ. However, when the GT operator 
commutes with the effective two-body interaction, the multiple commutators involved can be expressed as 
linear functions of the single particle energies allowing for an exact sum over all intermediate states [33 , 42, 43 . 
The Hamiltonian given by Eq. (13) does not include spin dependent two-body terms, allowing for the use of 
the summation method. While very convenient, and useful in the study of electromagnetic transitions [37], the 
model is explicitly excluding the spin-isospin interactions which are known to control the size of the 2u (3~ (3~ 
decay matrix elements in QRPA calculations [6] [18j [41]. Despite this strong limitation, PHFB calculations 
have proved to be able to reproduce the observed 2v (3~ (3~ decay half-lives through the explicit introduction of 
deformation [33]. An extended version of the present model, able to describe odd-odd nuclei, would be needed 
to include at the same time quadrupole-quadrupole and spin-isospin two body terms in the Hamiltonian, and 
evaluate the (3(3 decay matrix elements without employing the summation method. It represents a technical 
challenge which exceeds the objectives of the present work. Using the summation method, the NTME M 2v (k) 
can be written as 



M 2v (k) = — 
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where r„ is given by 



(7) 
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and 
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with 



£> A r M = [H, [H, 



(10) 



Assuming that the GT operator commutes with the effective two-body interaction, the Eq. p0[) can be further 
simplified to 

Eq + (k) 

where 

s(n„J„, j n ) - e(n v ,l v , j v ) for /3-/3"mode 



™ l 7 | e(»V, ^, jv) — e(i T , l-mi-n) for e + DBD modes v ' 

When a two-body isoscalar quadrupole-quadrupole interaction is employed, as in the case of the pseudo SU(3) 
model [39-41], the energy denominator is a well-defined quantity without any free parameter, because the GT 
operator commutes with this two-body interaction. In the pseudo SU(3) scheme the sum in the above equation 
reduces to a single term for 2v (3~ (3~ [42j|43] and 2v e + DBD modes [44) . 

In the present work, we use a Hamiltonian with PPQQ type [36] of effective two-body interaction. The 
Hamiltonian is explicitly written as 

H = H sp + V(P) + ( qq V(QQ) (13) 

where H sp denotes the single particle Hamiltonian. Further, ( qq denotes the strength of QQ part of the effective 
two-body interaction. The purpose of introducing it is to study the role of deformation by varying the strength 
parameter Q qq . The final results are obtained by setting the Q qq — 1. The pairing part of the effective two-body 
interaction V(P) is written as 

V{ p ) = - (?) TA-iy^-^-^aiaUsap (14) 



where a denotes the quantum numbers (nljm). The state a is same as a but with the sign of m reversed. The 
QQ part of the effective interaction V(QQ) is expressed as 



V(QQ) = -(f)£ £(- i rH^>0%-J ,5 > 44 a s a, (15) 



where 



16tt 



1/2 



fi=[^ L ) r*Y2(9,4>) (16) 



The model Hamiltonian given by Eq. (|13p is not isospin symmetric. Hence, the energy denominator is not as 
simple as in the case of pseudo SU(3) model. However, the violation of isospin symmetry for the QQ part of 
our model Hamiltonian is negligible, as will be evident from the parameters of the two-body interaction given 
later. Further, the violation of isospin symmetry for the pairing part of the two-body interaction is presumably 
small in the mass region under study. Under these assumptions, the expression to calculate the NTME My. v of 
2u e + DBD modes for + — > + transition in the PHFB model can be obtained by generalizing the above ideas, 
as follows. 

In the PHFB model, states with good angular momentum J are obtained from the axially symmetric HFB 
intrinsic state |$ ) with K=0 using the standard projection technique 45] given by 



l*oo> = 



(2J + 1) 



8tt 2 



D£ (n)R(n)\$ )<m (17) 



where J2(0) and Dq (Q) are the rotation operator and the rotation matrix respectively. The axially symmetric 
HFB intrinsic state |$o) can be written as 

\$ } = H(ui m + v im bl m bl fh )\0) (18) 

im 

where the creation operators b\ m and fot^ are defined as 

b\ m = £ ^«.™ a am an d ^\fa = £(~1)' + "' "^Cia.m^-m (19) 



The results of HFB calculations are summarized by the amplitudes (u; m , Wj m ) and expansion coefficients Cij >m , 
which are obtained through minimization of the expectation value of Hamiltonian given by Eq. (13) before 
angular momentum projection. While a Variation After Projection (VAP) has been shown to provide a more 
accurate description of the energy spectra |50j , they imply a numerical effort which seems to be unnecessary 
given the ability of the present version of the PHFB to reproduce the observed 2v (3~ f3~ decay matrix elements 

Finally, one obtains the following expression for the NTME M2 V (k) of the 2v j3 : /3 decay modes 



M 2v {k) = 
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where p and q in the last line refer to protons (neutrons) and neutrons (protons) in the parent (X) and daughter 
(Y) nuclei, for the (3~f3~ (e + DBD) mode, respectively. Further, 



n J = J {det[l + F^(9)f^}} 1 / 2 x {det[l + F^(9)f^]} l ' 2 d J Q0 {9) sm(6)d9 
o 



and 



n {m (0) = {det[l + F^(9)f^}}^ x {det[l + F^ (9)f^}}^ 
The matrices fy and Fx (0) are given by 

i 

and 
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The energy denominator is evaluated as follows. Within the present scheme the difference in proton and 
neutron single particle energies with the same quantum numbers explains the Isobaric Analog State energy, 
and is well described by the difference in Coulomb energies between the initial and intermediate nuclei. When 
appropriate, the spin-orbit energy splitting is also added. Hence 



^■{j^'Ki ^7r? Jit) ^ij^ut^ui ju) — 



A c for n u ^n^ : l u ^l^ : j lJ = j n 

A c + AE s . . spliuing for n u = n^,l u = ^ j n 

A c - 2E for n v =n. K ,l v = U,j v = j-w 
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where the Coulomb energy difference Ac is given by Bohr and Mottelson 

0.70 



A c = 



A 1 / 3 



(2Z + 1) - 0.76 Uz + 1) 4/3 - Z 4/3 } 



(25) 



(26) 



It must be underlined that, in the present context, the use of the summation method goes beyond the closure 
approximation, because each proton-neutron excitation is weighted depending on its spin-flip or non-spin-flip 
character. The explicit inclusion of the spin-orbit splitting in the energy denominator, Eq. (|25|) . implies that it 
cannot be factorized out of the sum in Eq. @. In this sense, employing the summation method in conjunction 



with the PHFB formalism is richer than what was done in previous application with the pseudo SU(3) model 

[alia]. 

To calculate the NTME M 2i , for the 2v [3f3 modes, the matrices fy and [Fx{8)] a p are evaluated using 
expressions given by Eqs. (]23[) and ([54"]) respectively. The required NTME M 2l/ (fc) are obtained using Eq. ([50"]) 
with 20 gaussian quadrature points in the range (0, 7r). 



3 Results and discussions 

The doubly even 100 Sn nucleus is treated as an inert core with the valence space spanned by 2s\ /2,l<^3/2> 1^5/2> 
I/7/2, O57/2, O/19/2 and O/in/2 orbits for protons and neutrons. The set of single particle energies (SPE's) 
used (in MeV) are e(ld 5/2 )=0.Q, e(2s 1/2 )=1.4 J e(ld 3/2 )=2.0, £(0# 7/2 )=4.0, e(0h 11/2 )=6.5 (4.8 for 150 Nd and 
150 Sm), e(l/ 7/2 )=12.0 (11.5 for 150 Nd and 150 Sm), e(0ft 9 / 2 )=12.5 (12.0 for 150 Nd and 150 Sm) for protons and 
neutrons. This set of SPE's, but for the e(0/i 11 / 2 ), which is increased by 1.5 MeV, has been employed in a 
number of successful variational model calculations for nuclear properties in the mass region ^4^130 [4"7] . The 
constant difference between proton and neutron SPE's with the same quantum numbers, due to the Coulomb 
interaction, is only relevant in the evaluation of the energy denominators, given by Eq. (25). The spin-orbit 
splitting AE S , , splitting is evaluated using the above mentioned SPE's. 

The strengths of the pairing interaction V(PP) is fixed through the relation G p = G n = 35/A MeV. The 
strengths of the like particle components of the QQ interaction are taken as: Xpp = Xnn = 0.0105 MeV b~ A , 
where b is the oscillator parameter. The strength of proton-neutron (pn) component of the QQ interaction Xpn 
is varied so as to obtain the yrast spectra in optimum agreement with the experimental data. The theoretical 
spectra is taken to be the optimum if the excitation energy of the 2 + state E 2 + is reproduced as closely as 
possible in comparison to the experimental results. Thus, we fix Xpn through the experimentally available energy 
spectra for a given model space, SPE's, G p , G n and Xpp- We present the values of Xpn in Table 1. These values 
of the strength of the QQ interaction are comparable to those suggested by Arima on the basis of an empirical 
analysis of the effective two-body interactions [IS] . All these input parameters are kept fixed throughout the 
calculation. Further, we have performed independent calculations for the parent and daughter nuclei involved 
in the j3(3 decay, whose deformations are in general different. 



3.1 Yrast spectra and electromagnetic properties 

In Table 1, we display the theoretically calculated and experimentally observed yrast spectra of 124 . 126 . 128 . 130 > 132 Xe, 
i24,i26,i 2 8,i30 Te i 30 ,i32 Ba i sotopes f or j« =2 + } 4+ an d 6+ states. All the experimentally observed E 2+ g§] 

energies are reproduced with two significant digits. However, it is noticed that the theoretical spectra is expanded 
in comparison to the experimental spectra. This effect can be corrected by complementing the PHFB model 
with the VAP prescription [50) . However, our aim is to reproduce properties of the low- lying 2 + states. Hence, 
we have not attempted to invoke the VAP prescription, which will unnecessarily complicate the calculations. 

We present the theoretically calculated as well as the experimentally observed reduced B(E2: + — > 2 + ) 
transition probabilities [51] . static quadrupole moments Q(2 + ) and the gyromagnetic factors g(2 + ) [52] in 
Table 2. In the case of B(E2: + — > 2 + ), only some representative experimental as well as adopted values 53J 
are tabulated. The calculated reduced B(E2: + — > 2 + ) transition probabilities are presented for three values 
of the effective charges, e e // =0.40, 0.50 and 0.60, in columns 2 to 4, respectively. The experimentally observed 
data are tabulated in column 5. It is observed that the calculated B(E2: + — > 2 + ) are in excellent agreement 
with the observed results in all cases except for 126 Xe and 150 Sm isotopes. The calculated and experimentally 
observed B(E2: + — » 2+) transition probability of 126 Xe and 150 Sm differ by 0.015 and 0.466 e 2 b 2 respectively 
for the same e e //. However, they are quite close for e e //=0.4. 

The theoretically calculated Q{2 + ) are displayed in columns 6 to 8 of Table 2, for the same effective charges 
as mentioned above. The experimental Q{2 + ) data [52] are given in column 9 of the same table. No experimental 
Q{2 + ) results are available for 124 < 12 M 2 8,i30x e anc i 132 B& nuclei. The calculated and experimental Q(2 + ) are 
in agreement for 150 Nd and 150 Sm nuclei. For rest of the nuclei, the theoretical Q(2 + ) results are quite off from 
the observed values. 

The g-factors g(2 + ) are calculated with gf =1.0, gf =0.0, and gl = g v s =0.60, and are presented in the 
column 10 of Table 2. The available experimental g(2 + ) data [55] are given in column 11. The calculated and 
experimentally observed g{2 + ) are in good agreement for 124 Te and 126,128 Xe isotopes, whereas they are off by 



0.052, 0.054, 0.065 and 0.013 nm only in case of 124 Xe, 126 Te, 130 Ba and 130 Xe isotopes respectively. However, 
the calculated g(2 + ) values of 128 < 130 Te, 132 Ba, 132 Xe, 150 Nd and 150 Sm nuclei are off from the experimental 
results. 

In general, the overall agreement between the calculated and observed electromagnetic properties is reason- 
ably good. From this analysis we conclude that the PHFB wave functions of i24,i26,i28,i30,i32 Xej i24,i26,i28,i30 Xei 
130 ' 132 Ba, 150 Nd and 150 Sm nuclei, generated by fixing Xpn to reproduce the energy of the first 2 + state, provide 
a reasonable input for calculating (3j3 decay half-lives. In the following, we present the results of the NTMEs 
M 2y , and the half-lives Tfl, of 128 < 130 Te, iso Nd nuclei for 2v p-p- decay, and 124 ' 126 Xe and 130 < 132 Ba nuclei 

for 2v c + DBD modes for the + — > + transition, using the same PHFB wave functions. 
3.2 Decay rates 

In the case of 2v P~ P~ decay, the phase space factors for the + -^0 + transition have been given by Doi et 
al. [3] and Boehm et al. for g A — 1-25 [5]. We take the phase phase space factors from Boehm et al. [5] for 
calculating half-lives of 2v P~ P~ decay of 128 . 130 Te and 150 Nd isotopes. We take the phase space factors for 
the 2v e + DBD modes from Doi et al. 15J when available. For the rest of 2v e + DBD emitters, we calculate 
the phase space factors following the prescription of Doi et al. [15], in the approximation C\ — 1.0, C 2 = 0.0, 
C3 = 0.0 and = R + i(e) + =1.0. The phase space integrals have been evaluated for g A = 1-25 

and 1.261, for the + — » + transition of 2v P~ P~ decay and 2v e + DBD modes, respectively [T5J. However, it 
is more justified to use the nuclear matter value of g A around 1.0 in heavy nuclei. Hence, the theoretical 2£% 
are calculated for both g A =\.§ and 1.261(1.25) for 2v e+DBD (2v P~~ P~ decay) modes. 

3.2.1 2v p-p- decay 

In Table 3, all the available experimental and theoretical results for 2v P~P~ decay of 128 > 130 Te and 150 Nd 
isotopes are presented. The 2v P~ P~ decay of 128 > 130 Te isotopes for the + — > + transition has been observed 
geochemically, and results from direct detection methods are also available in the case of 130 Te and 150 Nd 
isotopes. In the same Table 3, we also give the experimentally extracted values of NTMEs M 2v . In Table 4, we 
compile all the available experimental and theoretical results for 2v e + DBD modes of 124 ' 126 Xe, 130,132 Ba isotopes 
along with our calculated NTMEs M 2v , and the corresponding half- lives Ty 2 for the + — > + transition. 

In the case of 128 Te isotope, we take the phase space factor G 2l/ (P~ P~) = 8.475 xl0~ 22 y _1 . In the PHFB 
model, the calculated M 2v is close to the NTMEs extracted from experiments of Takaoka et al. [54] , Bernatovicz 
et al. [55J and Barabash [5S] for g A =1.0, while it is close to those of Lin et al. [57J and Manuel [55] for g A =1.25. 
The NTMEs M 2u calculated in SU(4) CTT [55], SSDH JBDJ and MCM gT] differ from M 2v calculated in the present 
work by factor of 1.4 - 1.6 and are in agreement with the extracted M 2v due to Lin et al. [57] and Manuel [55J 
for g A =1.0. The NTME M 2v calculated by Civitarese et al. [5J in SSDH is close to the NTMEs extracted from 
the experiments of Takaoka et al. [54], Manuel [63] and Barabash [56 for g A =1.25. The presently calculated 
NTME is smaller than the M 2v calculated in QRPA [64] and WCSM [3] by factor of approximately 2.2 and 
3.6 respectively while it is larger by a factor of approximately 5 than the M 2v calculated by Stoica [65] in 
SRPA(WS). The half-life T?Y 2 calculated in QRPA [H6] is close to the experimental Ty 2 due to Takaoka et al. 
[54] . Lin et al. [57] and Barabash [56]. The calculated half- lives JfY 2 by Hirsch et al. [67 in OEM and Scholten 
et al. [68] in IBM are quite small. In the present calculation, the predicted half-life T 2 ^ 2 is (1.05-2.55) x 10 24 y 
for g A = (1.25 - 1.0). 

The phase space factor G 2v {P~ ~/3~) used to study the 2v P~ /3~ decay of 130 Te is 4.808xl0" 18 y _1 . The 
NTME M 2v calculated in the PHFB model is in agreement with the NTME extracted from the experiment of 
Milano+INFN experiment 69J for g A =1.0. The NTME M 2v calculated by Caurier et al. [70] in the SM and 
Aunola et al. 6 lj in the MCM are close to the recommended value of Barabash [56] for g A =1.0. The NTME 
M 2v calculated by Stoica et al. 65J in SRPA (WS) is in agreement with the experimental M 2v of Lin et al. [57] . 
Takaoka et al. [54], Barabash [56] and Milano+INFN experiment 69J for g A =1.25. The presently calculated 
NTME differs from the NTME calculated in SU(4) CTr [5SJ and QRPA [53] by a factor of 1.1 and 1.2 respectively. 
In the RQRPA, the calculated M 2v by Toivanen et al. [71] is close to the experimental value of Bernatovicz et 
al. [55] for g A =1.25. The calculated half-lives Tff 2 by Hirsch et al. [67] in OEM, Scholten et al. [68] in IBM 
and Haxton et al. [3] in WCSM are quite small while the half-life Tfj 2 calculated by [TpJ in SU(4)o- T is close 



to the experimental results of Lin et al. [S7], Takaoka et al. [SI] and Milano+INFN experiment [55]. In the 
present calculation, the predicted half-life TfY 2 is (1.16-2.82) x 10 20 y for g A = (1.25 - 1.0). 

In the case of 150 Nd, the phase space factor G 2v {(3~~ (3~) used is 1.189xl0 -16 y _1 . The experimental M 2v 
of UCI [72] and NEMO 3 [73] for g A =1.25 and ITEP+INR [7J for g A =1.0 are close to the calculated value 
of M 2v in present work. The M 2v calculated by Hirsch et al. g3| in pseudo-SU(3) is in close agreement to the 
experimental M 2v of UCI 72J and NEMO 3 [73] for g A =1.0. The calculated value of M 2v by [42] in pseudo- 
SU(3) and [55] in SU(4) CTT is larger by a factor of 1.67 and 1.95 respectively than the presently calculated 
value of M 2v . The half-life T 2 ^ calculated by Hirsch et al. [57] in the OEM is close to the experimental result 
of ITEP+INR experiment while the Tffe calculated by Staudt et al. 66J in QRPA favors the value of UCI 
experiment. The predicted half-life T\% in our PHFB model is (7.89-19.27) x 10 18 y for g A = (1.25 - 1.0). 

3.2.2 2v e+DBD modes 

In the case of 124 Xe, the available experimental half-life limits on Tfj 2 of 2v f3 + (3 + and 2v EC EC modes, for the 
+ — » + transition, are > 8.0 x 10 23 y [84] and > 1.1 x 10 17 y [85] respectively. Barabash [86] has investigated 
2v (3 + (3 + and 2v (3+EC modes and the half-life limits T 2 ^ are > 2.0 x 10 14 y and > 4.8 x 10 16 y respectively. 
The phase space factors used in the present calculation are G 2v {(3+ (3+) = 1.205xl0 -25 y _1 , G 2v ((3+EC) = 
4.353xl0~ 21 y- 1 and G 2v {ECEC) = 5.101xl0~ 20 y" 1 . The calculated NTMEs M 2u in the PHFB and SU(4) CTT 
[59] models are close to each other for 2v (3+ EC and 2v ECEC modes. The presently calculated NTME M 2v 
is larger than the NTME due to Aunola et al. [5T] in MCM by a factor of 5 and 7 approximately for 2v (3+EC 
and 2v ECEC modes respectively while the NTMEs M 2v calculated in the QRPA model [57] are larger than 
that of PHFB model values by a factor of 1.5 approximately for all the three modes. The calculated half- life 
Ty 2 by Staudt et al. [88 differ by two orders of magnitude form the half-lives Tfj 2 calculated in the present 
work for all modes. The theoretically calculated Ty 2 are of the order of 10 25-27 y, io 22-24 y and 10 21-24 y for 
2v (3+/3+, 2v (3+EC and 2v ECEC modes respectively for g A = 1.261 - 1.00. The calculated half-lives Ty 2 
in the PHFB model are (3.015-7.624) x 10 27 y, (8.347-21.11) x 10 22 y and (7.123-18.01) x 10 21 y for 2v (3+ (3+ , 2v 
/3+EC and 2v ECEC modes respectively for g A =(1.26-1.0). 

The e + DBD modes of 126 Xe isotope for the + — » + transition have not been studied so far either experi- 
mentally or theoretically. In the present calculation, the predicted half-life Ty 2 is (5.682-14.37) xlO 24 y for g A 
= (1.261-1.0) using the phase space factor G 2v {ECEC) = 7.428x 10~ 23 y" 1 for 2v ECEC mode. 

In the case of 130 Ba isotope, the + — » + transition of c + DBD modes has been investigated in geochemical 
experiments by Meshik [55] as well as Barabash [50J and combined limits for all the e + DBD modes are given. 
Meshik [S5J has given a half-life limit Tfj 2 = (2.16 ± 0.52) x 10 21 y. The half-life limit for (2v + Ou) decay 

mode given by Barabash [90] is > 4.0 x 10 21 y for all the c + DBD modes. The phase space factors used in the 
present calculation are G 2u {(3+ P+) = 1.211xl0- 27 y" 1 , G 2v (/3 + EC) = 1.387xl0- 21 y" 1 and G 2v {ECEC) = 
4.134xl0~ 20 y- 1 respectively. The theoretical NTME M 2lJ calculated in the PHFB and SU(4) CTT [S5] models 
differ by a factor of -1.4 for 2v /3+EC and 2v ECEC modes. The NTME M 2v calculated in the PHFB model is 
smaller than the NTMEs of Hirsch et al. [57] by a factor of 1.7 in the case of 2v (3 + [3 + mode while for 2v (3 + EC 
and 2v ECEC modes, the results differ by a factor of 2 approximately. On the other hand, the calculated 
NTME in the present work is larger than the calculated M 2v due to Aunola et al. [5T] in MCM by a factor 
of 2.7 and 6 approximately for 2v (3 + EC and 2v ECEC modes respectively. The calculated half-life T 2 ^ 2 by 
Staudt et al. ;SS] in QRPA model for 2v /3 + /3 + mode is of the order of present calculation. The predicted Ty 2 
oi2v f3 + f3+, 2v (3+EC and 2v ECEC modes in PHFB model are (4.797 - 12.13) x 10 29 y, (4.188 - 10.59) x 10 23 
y and (1.405 — 3.553) xlO 22 y respectively for g A — (1.26 — 1.0). The experimentally observed half-life for the 
0+ ^0+ transition of 130 Ba [55] is close to the theoretically predicted Tfj 2 of 2v ECEC mode. 

For the 132 Ba isotope, the deduced total half-life limit for of all the e + DBD modes (2v + Qv decay) for the 
+ — » + transition by Barabash [50] in geochemical experiment is > 3.0 x 10 20 y. Meshik [55J has also carried 
out the geochemical experiment and has given a combined half-life limit Ty 2 — (1.3 ± 0.9) x 10 21 y for all the 
e + DBD modes. No theoretical study has been done so far for the study of e + DBD modes of 132 Ba. The phase 
space factor G 2w (ECEC) used in the present calculation is 6.706xl0 -23 y _1 . In the present calculation, the 
predicted half-life T 2l/ 2 for the 2v ECEC mode is (5.474 - 13.84) x 10 24 y for g A = (1.261 - 1.0). 



3.3 Deformation effect 



There are several possibilities to quantify the deformation parameter of the nucleus. We take the quadrupole 
moment of the intrinsic state (Qq) (in arbitrary units) and the quadrupole deformation parameter 02 as a 
quantitative measure of the deformation. The variation of the (Qq), P2 and M.2 V with respect to the change 
in the strength of the QQ interaction ( qq has been investigated to understand the role of deformation on the 
NTME M 2 v In table 5, the quadrupole moment of the intrinsic states (Qq), deformation parameter 02 and the 
NTMEs M.2 V for different ( qq are presented. We calculate the deformation parameter with the same effective 
charges as used in the calculation of reduced B(E2:0 + — > 2 + ) transition probabilities. 

It is noticed that the (Q 2 ,) as we ^ as increases in general as the ( qq is varied from to 1.5 except a few 
anomalies. In all cases, it is found that the quadrupole deformation parameter /3 2 follows the same behavior 
as the quadrupole moment of the intrinsic state (Qq) with respect to the change in ( qq . Further, there is an 
anticorrelation between the deformation parameter and the NTME Af 2 „ in general but for a few exceptions. 

We define a quantity _D 2t , as the ratio of Mi v at zero deformation = 0) and full deformation (C, qq = 1) 
to quantify the effect of deformation on M 2 j,. The ratio Z? 2i , is given by 

E>2u - 7T 7T ■ ( Z 

M 2 v{Q q = 1) 

The ratios D 2u are 3.63, 3.47, 4.26, 2.89, 4.66, 3.10 and 6.02 for 124 Xe, 126 Xe, 128 Te, 130 Te, 130 Ba, 132 Ba 
and 150 Nd nuclei respectively. These values of D 2 j, suggest that the M 2v is quenched by a factor of 3 to 6 
approximately in the mass region A ~ 120 — 150 due to deformation effects. In contrast to other models, the 
quenching of the NTMEs seems to be closely related with the explicit inclusion of deformation effects. 

In (3/3 decay studies where deformation was included but no angular momentum projection was performed, 
nuclear deformation was found to be a mechanism of suppression of the 2is (3(3 decay. In this case, the (3(3 decay 
matrix elements are found to have maximum values for about equal deformations of parent and daughter nuclei, 
and they decrease rapidly when differences in deformations increase [32]. This deformation effect is different 
from the one reported in this work. Further research is needed to relate these two approaches. 

The quenching factors discussed above could be considered as a conservative estimate of the uncertainties 
in the predicted NTMEs of 2v (3(3 modes, given the fact that many of the nuclei studied are in the transitional 
region and do not display a well defined rotational spectrum. Also the PPQQ interaction employed in the 
present calculation is of schematic nature. These uncertainties qualify both the present results, those obtained 
with other models where deformation is not explicitly considered, or where the rotational symmetry is not 
restored. The uncertainties associated with the 0v (3(3 decay processes would be expected to be far smaller than 
in the 2v (3(3 modes. 



4 Conclusions 

To summarize, we have built the PHFB wave functions of the parent and daughter nuclei involved in (3(3 decay 
processes. Theoretical results for the yrast spectra, reduced B(E2: + — * 2 + ) transition probabilities, static 
quadrupole moments Q(2+) and ^-factors g{2+) of i24,i26,i28,i30,i32 Xe; 124,126,128,130^ 130,132^ i 5 o Nd and 

150 Sm nuclei were presented and compared with the available experimental results. The same PHFB wave 
functions were employed to calculate NTMEs M 2u and half-lives T\ v /2 of 124 Xe (2v (3+ (3+ , 2v (3+EC and 2v 
ECEC modes), 126 Xe {2v ECEC mode), 128 Te {2v (3~ (3~ decay), 130 Te {2v (3~ (3~ decay), 130 Ba {2v (3+ (3+ , 
2v (3+ EC and 2v ECEC modes), 132 Ba {2v ECEC mode) and 150 Nd (2v (3' (3~ decay) isotopes. It was shown 
that the proton-neutron part of the PPQQ interaction, which is responsible for triggering deformation in the 
intrinsic ground state, plays an important role in the quenching of Mzv, by a factor of approximately 3 to 6, in 
the considered mass region A ^120-150. 

This work was partially supported by DAE-BRNS, India vide sanction No. 2003/37/14/BRNS/669, by 
Conacyt-Mexico and DGAPA-UNAM. 



References 

[1] H. Primakoff, Phys. Rev. 85, 888 (1952); S.P. Rosen, H. Primakoff, Report Prog, in Phys. 22, 121 (1959); 
S.P. Rosen, H .Primakoff, Alpha-beta-gamma ray spectroscopy, ed. K. Siegbahn (1965) p. 1499; H. Primakoff, 
S.P. Rosen, Ann. Rev. Nucl. Part. Sci. 31, 145 (1981). 



[2 
[3 

[4; 

[5 

[6 

[7: 

[8 

[9 
[10 

[11 
[12 

[13 

[14; 

[15 
[16 
[17 
[18 
[19 
[20 
[21 

[22 
[23 
[24 
[25 
[26 

[27 
[28 



D. Bryman, C. Picciotto, Rev. Mod. Phys. 50, 11 (1978). 

W. C. Haxton, G. J. Stephenson, Jr., Prog. Part. Nucl. Phys. 12, 409 (1984). 
M. Doi, T. Kotani, E. Takasugi, Prog. Theo. Phys. Suppl. 83, 1 (1985). 

J. D. Vergados, Nucl. Phys. B 218, 109 (1983), Phys. Rep. 133, 1 (1986); ibid 361, 1 (2002). 
A. Faessler, Prog. Part. Nucl. Phys. 21, 183 (1988). 
T. Tomoda, Rep. Prog. Phys. 54, 53 (1991). 

F. Bochm, P. Vogel, Ann. Rev. Nucl. Part. Sci. 34, 125 (1984); F. Boehm, P. Vogel, Physics of Massive 
Neutrinos, 2nd ed. (Cambridge University Press, Cambridge 1992). 

K. Zuber, Phys. Rep. 305, 295 (1998). 

E. Fiorini, Phys. Rep. 307, 309 (1998). 
H. Ejiri, Phys. Rep. 338, 265 (2000). 

S. R. Elliott and J. Engel, J. Phys. G: Nucl. Part. Phys. 30 (2004) R183. 

V.I. Tretyak, Y.G. Zdesenko, At. Data Nucl. Data Tables 61, 43 (1995); ibid 80, 83 (2002). 

H.V. Klapdor-Kleingrothaus, et al., Mod. Phys. Lett. 16 (2001) 2409; H.V. Klapdor-Kleingrothaus, et al., 
Phys. Lett. B 586 (2004) 198. 

M. Doi, T. Kotani, Prog. Theor. Phy. 87, 1207 (1992). 

M. Doi, T. Kotani, Prog. Theor. Phy. 89, 139 (1993). 

M.K. Moe, P. Vogel, Ann. Rev. Nucl. Part. Sci. 44, 247 (1994). 

J. Suhonen, O. Civitarese, Phys. Rep. 300, 123 (1998). 



A. Faessler, F. Simkovic, hep-ph/9901215; J. Phys. G 24, 2139 (1998). 

H.V. Klapdor-Kleingrothaus, Sixty years of Double Beta Decay, World Scientific, Singapore (2001). 

A.S. Barabash, V.I. Umatov, R. Gurriaran, F. Hubert, Ph. Hubert, M. Aunola, J. Suhonen, Nucl. Phys. 
A 604, 115 (1996). 

O. Civitarese, J. Suhonen, Phys. Lett. B 482, 368 (2000). 
J. Suhonen, O. Civitarese, Phys. Lett. B 497, 221 (2001). 

A. Griffiths, P. Vogel, Phys. Rev. C 46, 181 (1992). 

J. Suhonen, O. Civitarese, Phys. Rec. C 49, 3055 (1994). 

N. Auerbach, D.C. Zheng, L. Zamick, B.A. Brown, Phys. Lett.B 304 (1993) 17; N. Auerbach, G.F. Bertsch, 

B. A. Brown, L. Zhao, Nucl. Phys. A 56, 190 (1993). 

D. Troltcnicr, J.P. Draayer, J.G. Hirsch, Nucl. Phys. A 601, 89 (1996). 
F. Frisk, I. Hamamoto, X.Z. Zhang, Phys. Rev. C 52, 2468 (1995). 



[29] P. Sarriguren, E. Moya de Guerra, A. Escuderos, A.C. Carrizo, Nucl. Phys. A 635, 55 (1998); P. Sarriguren, 
E. Moya de Guerra, A. Escuderos, Nucl. Phys. A 658 (1999) 13; ibid Nucl. Phys. A 691, 631 (2001). 



[30 
[31 
[32 

[33 
[34 

[35 
[3G 
[37" 
[38 
[39 

[40 
[41 
[42 
[43 
[44" 
[45 
[46 
[47" 
[48 
[49 
[50 

[51 

[52 
[53 
[54; 
[55 



E. Nacher et al., Phys. Rev. Lett. 92, 232501(2004). 

L. Pacearescu, A. Faessler, F. Simkovic, Phys. At. Nucl. 67, 1210 (2004). 

R. Alvarez-Rodriguez, P. Sarriguren, E. Moya de Guerra, L. Pacearescu, A. Faessler, F. Simkovic, Phys. 
Rev. C 70, 064309 (2004). 

R. Chandra, J. Singh, P.K. Rath, P.K. Raina, J.G. Hirsch, Eur. Phys. J. A 23, 223 (2005) 

P.K. Raina, A. Shukla, S. Singh, P.K. Rath, J.G. Hirsch, Eur. Phys. J. A 28, 27 (2006); A. Shukla, P.K. 
Raina, R. Chandra, P.K. Rath, J.G. Hirsch, Eur. Phys. J. A 23, 235 (2005). 

A. Shukla, P.K. Raina, P.K. Rath, J.Phys. G: Nucl. Part. Phys. 34, 549 (2007). 
M. Baranger, K. Kumar, Nucl. Phys. A 110, 490 (1968). 

B. M. Dixit, P.K. Rath, P.K. Raina, Phys. Rev. C 65, 034311 (2002), Phys. Rev. C 67, 059901(E) (2003). 
O. Civitarese, J. Suhonen, Phys. Rev. C 47, 2410 (1993). 

X. R. Wu, A. Staudt, H. V. Klapdor-Kleingrothaus, Chen-rui Ching and Tso-hsiu Ho, Phys. Lett. B 272, 
169 (1991); X. R. Wu, A. Staudt, T. T. S. Kuo and H. V. Klapdor-Kleingrothaus, Phys. Lett. B 276, 274 
(1992). 

J. Engel, W. C. Haxton and P. Vogcl, Phys. Rev. C 46, R2153 (1992) 
P. Vogel, M.R. Zirnbauer, Phys. Rev. Lett. 57, 3148 (1986). 

O. Castahos, J.G. Hirsch, O. Civitarese, P.O. Hess, Nucl. Phys. A 571, 276 (1994). 
J.G. Hirsch, O. Castahos, P.O. Hess, O. Civitarese, Phys. Rev. C 51, 2252 (1995). 
V.E. Ccron, J.G. Hirsch, Phys. Lett.B 471, 1 (1999). 
N. Onishi, S. Yoshida, Nucl. Phys. A 260, 226 (1966). 

A. Bohr, B.R. Mottelson, Nuclear Structure Vol. I (World Scientific, Singapore, 1998). 
Rani Devi, S.P. Sarswat, Arun Bharti, S.K. Khosa, Phys. Rev. C 55, 2433 (1997). 
A. Arima, Nucl. Phys. A 354, 19 (1981). 
M. Sakai At. Data, Nucl. Data Tables 31 , 400(1984). 

S.K. Khosa, P.N. Tripathi, S.K. Sharma, Phys. Lett. B 119, 257 (1982); P.N. Tripathi, S.K. Sharma, S.K. 
Khosa, Phys. Rev. C 29, 1951 (1984); S.K. Sharma, P.N. Tripathi, S. K. Khosa, Phys. Rev. C 38, 2935 
(1988). 

S. Raman, C.H. Malarkey, W.T. Milner, CW. Nestor, Jr., P.H. Stelson, At. Data Nucl. Data Tables 36, 1 
(1987). 

P. Raghavan, At. Data Nucl. Data Table 42, 189 (1989). 

S. Raman, C. W. Nestor, Jr., P. Tikkanen, At. Data Nucl. Data Tables 78 (2001). 
N. Takaoka et al, Phys. Rev. C 53, 1557 (1996). 

T. Bernatovicz, J. Brannon, R. Brazzle, R. Cowsik, C. Hohenberg, F. Podosek, Phys. Rev. Lett. 69, 2341 
(1992); Phys. Rev. C 47, 806 (1993). 



[56] A.S. Barabash, Czech. J. Phys. 52, 567 (2002); nucl-ex/0203001 (2002). 



[57] J.W. Lin et ai, Nucl. Phys. A 481, 484 (1988). 

[58] O.K. Manuel, Proc. Int. Symp. on Nuclear Beta Decays, Neutrino, Osaka, Japan, June 1986 (World Sci., 
1986), pg. 71. 

[59] O.A. Rumyantsev, M.H. Urin, Phys. Lett. B 443, 51 (1998). 

[60] S.V. Semenov, F. Simkovic, V.V. Khruschev, P. Domin, Phys. At. Nucl. 63, 1196 (2000). 
[61] M. Aunola, J. Suhonen, Nucl. Phys. A 602, 133 (1996). 
[62] O. Civitarese, J. Suhonen, Phys. Rev. C 58, 1535 (1998). 
[63] O.K. Manuel, J. Phys. G17, S221 (1991). 

[64] J. Engel, P. Vogel, M.R. Zirnbauer, Phys. Rev. C 37, 731 (1988). 
[65] S. Stoica, Phys. Rev. C 49, 2240 (1994). 

[66] A. Staudt, K. Muto, H.V. Klapdor, Euro. Phys. Lett. 13, 31 (1990). 

[67] M. Hirsch, X.R. Wu, H.V. Klapdor-Klcingrothaus, Ching Cheng-rui, Ho Tso-hsiu, Phys. Rep. 242, 403 
(1994). 

[68] O. Scholten, Z.R. Yu, Phys. Lett. B 161, 13 (1985). 
[69] C. Arnaboldi et at, Phys. Lett. B 557, 167 (2003) 

[70] E. Caurier, F. Nowacki, A. Poves, J. Retamosa, Nucl. Phys. A 654, 973 (1999). 
[71] J. Toivancn, J. Suhonen, Phys. Rev. C 55, 2314 (1997). 

[72] A. De Silva, M.K. Moe, M.A. Nelson, M.A. Vient, Phys. Rev. C 56, 2451 (1997). 

[73] D. Lalanne, [hep-ex/0509005] Dans High Energy Physics ICHEP 2004 (in 2 volumes), World Scientific (Ed.) 
(2005). 

[74] V. Artemiev et ai, Phys. Lett. B 345, 564 (1995). 

[75] T. Kirsten et ai, Proc. Int. Symp. "Nuclear Beta Decay, Neutrino (Osaka 86)", World Sci., Singapore, 81 
(1986). 

[76] A. Allessandrello et ai, Phys. Lett. B 486, 13 (2000). 

[77] E. Bellotti, C. Cattadori, O. Cremonesi, E. Fiorini, C. Liguori, A. Pullia, P.P. Sverzellati, L. Zanotti, Euro. 
Phys. Lett. 3, 889 (1987). 

[78] Yu. G. Zdesenko, Sov. J. Part. Nucl. 11(6), 542 (1980). 

[79] O.A. Rumyantsev, M.H. Urin, JETP Lett. 61, 361 (1995). 

[80] S.I. Vasilev et ai, JETP Lett. 51, 622 (1990); 58, 178 (1993). 

[81] V. Artemiev et ai, JETP Lett. 58, 262 (1993). 

[82] S.R. Elliott, M.K. Moe, M.A. Nelson, M.A. Vient, Nucl. Phys. B (Proc. Suppl.) 31, 68 (1993). 

[83] A. A. Klimenko, A. A. Pomansky, A. A. Smolnikov, Nucl. Instrum. Meth. B 17, 445 (1986). 

[84] A. Bolozdynya et. ai, IEEE Nuclear Science Symposium, Medical Imaging Conference, 2, 697 (1996); ibid 
IEEE Transaction on Nuclear Science, 44, 3 (1997). 

[85] Ju.M. Gavriljuk, V.V. Kuzminov, NYu. Osetrova, S.S. Ratkevich, G.V. Volchenko, Phys. At. Nucl. 61, 
1287 (1998). 



[86] A.S. Barabash, V.V. Kuzminov, V.M. Lobashev, V.M. Novikov, B.M. Ovchinnikov, A. A. Pomansky, Phys. 
Lett. B 223, 273 (1989). 

[87] M. Hirsch, K. Muto, T. Oda, H.V. Klapdor- Klcingrothaus, Z. Phys. A 347, 151 (1994). 

[88] A. Staudt, K. Muto, H.V. Klapdor-Kleingrothaus, Phys. Lett. B 268, 312 (1991). 

[89] A.P. Mcshik, CM. Hohcnbcrg, O.V. Pravdivtseva, Ya. S. Kapusta, Phys. Rev. C 64, 035205 (2001). 

[90] A. S. Barabash, R.R. Saakyan, Phys. At. Nucl. 59, 179 (1996). 



Table 1: Excitation energies (MeV) of J n = 2+ 4+ and 6+ yrast states of 124,126,128,130^ i 24 ,i26,i28,i30,i32 Xe) 
and 130 ' 132 Ba, 150 Nd and 150 Sm nuclei. 
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0.4645 


0.4646 




0.01536 


E 2 + 


0.6675 


0.6677 






E 4 + 


1.4519 


1.1277 






E 4+ 


1.6278 


1.4403 






E 6 + 


2.7982 


1.9328 






E 6 + 


2.7263 


2.1118 


150 Nd 


0.02160 


E 2 + 


0.1307 


0.13012 


15U Sm 


0.01745 


E 2 + 


0.3328 


0.33395 






E 4 + 


0.4320 


0.3815 






E 4+ 


1.0156 


0.77335 






E 6 + 


0.8960 








E 6 + 


1.9185 


1.27885 



Table 2: Comparison of theoretically calculated and experimentally observed reduced transition probabilities 
B(E2: 0+ -> 2+), static quadrupole moments Q(2+) and .9 factors g( 2+) f 124 > 126 . 128 < 130 Te, i24,i26,i28,i30,i32 Xej 
and 130 > 132 Ba nuclei. Here B(E2) and Q{2+) are calculated for effective charge e^ — l-\-e e ff and e n — &eff- * 
denotes the average B(E2) values. 



Nucleus 




B(E2 : 0+ 2+) (e 2 b 2 ) 
Theory Exp. [ST] 




Q(2+) (cb) 
Theory 


Exp.[SJ] 


5(2 
Theory 


"*") (nm) 
Exp.JSJ 






e e// 








e e// 












0.40 


0.50 


0.60 




0.40 


0.50 


0.60 








124 Xe 


0.722 


0.941 


1.188 


qfi+n n«* 

. <^vj _ i— \J . \J\J 

0.90±0.07 
1.49±0.09 


-0.770 


-0.879 


-0.988 




0.302 


0.23±0.02 


l24 Te 


0.401 


0.530 


0.677 


5fis+n no6* 

0.61±0.20 
0.538±0.028 


-0.574 


-0.659 


-0.745 


-0.45±0.05 


0.348 


0.33±0.03 
0.31±0.04 
0.28±0.03 




0.670 


0.865 


1.084 


n 770+n f)25* 

0.79±0.06 

0.760±0.026 


-0.742 


-0.843 


-0.943 




0.373 


0.37±0.07 
0.27±0.04 




0.353 


0.460 


581 


0.487±0.035 
0.532±0.037 


-0.539 


-0.614 


-0 691 


-0.20±0.09 


0.424 


0.34±0.03 
0.31±0.04 


l2 «Te 


0.298 


0.381 


474 


^«^+n oofi* 

yj .UOO- 1— \J. \J\J\J 

0.380±0.009 
0.387±0.011 


-0.496 


-0.561 


-0 696 


-0.06±0.05 
-0.14±0.12 


0.514 


0.35±0.04 
0.31±0.04 


l2 «Xe 


0.637 


0.819 


1.024 


7c,n+o 040* 

0.790±0.040 
0.890±0.230 


-0.724 


-0.820 


-0 91 7 




0.400 


0.41±0.07 
0.31±0.03 


iau Te 


0.231 


0.289 


0.354 


295±0 007* 

0.290±0.011 

0.260±0.050 


-0.438 


-0.490 


-0.542 


-0.15±0.10 


0.679 


0.33±0.08 
0.29±0.06 


13,, Xe 


0.493 


0.624 


0.769 


0.65±0.05* 

0.631±0.048 

0.640±0.160 


-0.637 


-0.716 


-0.795 




0.463 


0.38±0.07 
0.31±0.04 


13lJ Ba 


1.048 


1.331 


1.649 


1.163±0.016* 

1.36±0.14 

1.21±0.38 


-0.928 


-1.046 


-1.164 


-0.86±0.08 
-0.33±0.24 


0.445 


0.35±0.03 


132 Ba 


0.624 


0.767 


0.926 


0.86±0.06* 

0.73±0.18 

0.86±0.06 


-0.718 


-0.796 


-0.874 




0.649 


0.34±0.03 


132 Xe 


0.329 


0.409 


0.498 


0.460±0.030* 

0.42±0.11 

0.44±0.03 


-0.519 


-0.579 


-0.639 


0.010±0.005 


0.575 


0.39±0.05 
0.37±0.05 


150 Nd 


2.132 


2.580 


3.070 


2.760±0.040* 

2.640±0.080 

2.670±0.100 


-1.322 


-1.455 


-1.587 


-2.00±0.51 


0.636 


0.422±0.039 
0.322±0.009 



Sm 1.707 2.056 2.437 1.350±0.030* -1.182 -1.297 -1.412 -1.32±0.19 0.592 0.385±0.027 

1.470±0.090 -1.25±0.20 0.411±0.032 

1.440±0.150 



*data taken from Reference [53] 



Table 3: Experimental half- lives T^ 2 and corresponding matrix element M%v of 1v (3 /3 decay for the 
+ — > + transition of 128 > 130 Te and 150 Nd nuclei along with the theoretically calculated Mi v in different 



models. The numbers corresponding 


to (a) and (b) 


are 


calculated for 


A = 


1.25 and 1.0 r 


espectively. 






Nuclei 




Experiment 










Theory 








Ref. 


Project 


rp2v 

1/2 




\M 2v \ 


Ref. 


Model 


\M 2v \ 




rp2v 

1/2 


i28 Te 


[541 


gch. 


(2.2±0.3) 




023 + „ 00 ,°,^ 


* 


PHFB 


0.033 


fal 


1.05 


(10 24 y) 








(b) 


036 + n ^ 

u.uou_ 0022 








fbl 


2.55 




m 


gch. 


2.0 


(el) 
(b) 


0.024 
0.038 


60 


SSDH 


0.048 


fal 
(b) 


0.51 
1.29 




m 


Ech 

o 


7.7±0.4 


(a) 


n m n+0.0003 
U.U1Z_ 0003 


123 


SM 






0.5 










(b) 


n ni q+ 0005 

u - uly -0.0005 




SU(4) CTr 


0.053 


(a) 


0.42 




m\ 


gch. 


1.8±0.7 


(a) 


026+ 0071 

U.UZU_ 0039 








(b) 


1.06 










(h) 


040+ 0111 

U.U1U_ 0060 




SSDH 


013 


I'm) 

(A) 


6 98 




m 


gch. 


>5.0 


(a) 


<0.015 








(b) 


17.65 










(b) 


<0.024 




MCM 


0.046 


(a) 


0.56 




[58] 


gch. 


1.4±0.4 


(a) 


029+ 0053 

U - UZJ -0.0034 








(b) 


1.41 










(b) 


U.U4O_ 0053 


[65] 


SRPA(WS) 


0.006 


(a) 


32.78 




m 


Average 


7.2±0.3 


(a) 


n m q+o.0003 

U.UlO_ 0003 








(b) 


82.87 






Value 




(b) 


020+ 0004 

U.UZU_ 0004 


[67] 


OEM 






0.21 




m 


Recommended 


2.5±0.4 


(a) 


n n99 +0.0020 
u.uzz_ 0015 


[55] 


QRPA 






2.63 






Value 




(b) 


34+ 0031 

U.UM_ 0024 


[sa 

EH] 
[3] 


QRPA 

IBM 

WCSM 


0.074 
0.120 


(a) 
(b) 

(a) 
(b) 


0.22 
0.54 
0.09 
0.08 
0.21 



130 Te 
(10 20 yr) 





Milano+INFN 


6.1±1.4± 2 : 9 5 


(a) 


n ni q+o.0232 

U.UlS_ .oo43 


* 


PHFB 


0.042 


(a) 


1.16 






(b) 


029+ 0362 

u.uzy_ 0068 








(b) 


2.82 


m 


Milano 


>3.0 


(a) 
(b) 


<0.026 
<0.041 




SM 


0.030 


(a) 
(b) 


2.3 
5.84 


[54] 


gch. 


7.9±1.0 


(a) 
(b) 


016+ 0011 

U.U1D_ 0009 

025+ - 0018 
u.uzo_ 0015 


[59] 


SU(4) CTr 


0.0468 


(a) 
(b) 


0.95 
2.32 


[55] 


gch. 


27.0±1.0 


(a) 
(b) 


09 + 0002 
u.uua_ 0002 

n 01 4+0 0003 

U.U11_ 0002 


m 


RQRPA(AWS) 


0.009 


(a) 
(b) 


25.68 
62.70 


[63] 


gch. 


8.0 


(a) 
(b) 


0.016 
0.025 


m 


RQRPA(WS ) 


0.009 


(a) 
(b) 


25.68 
62.70 


m\ 


gch. 


7.5±0.3 


(a) 
(b) 


n m 7+0.0003 

' -0.0003 

026+ 0005 
u.uzd_ 0005 


m 


MCM 


0.028 


(a) 
(b) 


2.65 
6.48 


m\ 




>8.0 


(a) 


< 0.016 




SU(4) CTr 






7.0 








(b) 


<0.025 


m 


SRPA(WS) 


0.016 


(a) 


8.12 


m 


INR Kiev 


>0.0001 


(a) 
(b) 


<4.563 
<7.130 




OEM 




(b) 


19.84 
0.79 


m 


Average 


27±1.0 


(a) 


009+ 0002 

U.UU»_ 0002 


66J 


QRPA 






18.4 




Value 




(b) 


oi4+o o° 03 

U.U11_ 0002 


m 


QRPA 


0.049 


(a) 


0.87 


m 


Recommended 


9.0±1.5 


(a) 


01 5+o o° 15 
u.uio_ ooil 








(b) 


2.12 




Value 




(b) 


024+ - 0023 

U.UZ1_ 0018 


[68] 
® 


IBM 

WCSM 


0.114 


(a) 
(b) 


0.17 
0.16 
0.40 



Table 3 continued 



02 o+().()o3o 

u - UZM -0 .0023 

046 + 0047 
u.u^u_ 0036 

U.UOO_q 0025 
+.00051 

u.uoo_ 0038 
021 + 0043 

U.UZ1_ Q 003 g 
n oo+0.0067 

U.UOO_Q QQ 5 g 

<0.028 
<0.043 
n n99+ 0092 

u - uzz -0.0O56 

o os c i +0 - 0144 

0.031 

0.048 

<0.022 

<0.034 

u.uoo_ 0135 

u.uo4_ 0211 

n n^+ - 0052 
u.uoo_ 0036 

U ' UO4 -0.0056 



150 Nd 

(10 18 y) 



[73] NEMO 3 
[72] UCI 



9.7±0.7±1.0 



6.75t°il ± 0.68 



[74] ITEP +INR 18.8±|;| ± 1.9 



>11.0 



[S] ITEP +INR 17l5 U ±3.5 
9.0 



[82] UCI 
[33] INR 

[12] Average 
Value 

[56] Average 
Value 



>18 



7 0+ 11 ' 8 

'■ u -0.3 



7.0±1.7 



(a) 
(b) 
(a) 
(b) 
(a) 
(b) 
(a) 
(b) 
(a) 
(b) 
(a) 
(b) 
(a) 
(b) 
(a) 
(b) 
(a) 
(b) 



* 


PHFB 


0.033 


(a) 


7.89 








(b) 


19.27 


m 


SU(4) ffr 


0.0642 


(a) 


2.04 








(b) 


4.98 


m 


P SU(3) 


0.055 


(a) 


2.78 








(b) 


6.79 


m 


OEM 






16.6 


66J 


QRPA 






7.37 



gch. denotes the geochemical experiment, * denotes the present work 



Table 4: Experimental limits on half-lives Tfj 2 (0+ — > + ), theoretically calculated M2„ and corresponding 
7^(0+ -» 0+) for 2^ c+DBD modes of 124 < 126 Xe and 130 < 132 Ba nuclei. The numbers corresponding to (a) and 
(b) are calculated for ^=1.261 and 1.0 respectively. 



Nuclei 


Decay 




Experiment 






Theory 








Mode 


Ref 


^ (y) 


Ref. 


Model 


\M 2u \ 




t% 2 (y) 


i24 Xe 


f3+(3+ 


m 
m 


>8.0xl0 23 
>2.0xl0 14 


* 


PHFB 


0.0525 


(a) 
(b) 


3.015xl0 2Y 
7.624xl0 27 










[87] 


QRPA 
QRPA 


0.0770 


(a) 
(b) 


1.400xl0 27 
3.539xl0 27 
8.170xl0 25 




(3+ EC 


[86] 


>4.8xl0 lb 


* 


PHFB 


0.0525 


(a) 
(b) 


8.347xl0 22 
2.111xl0 23 












SU(4) CTT 


0.0528 


(a) 
(b) 


8.240xl0 22 
2.084xl0 23 












MCM 


0.0100 


(a) 
(b) 


2.297xl0 24 
5.809xl0 24 










[57] 


QRPA 


0.0875 


(a) 
(b) 


3.001 xlO 22 
7.587xl0 22 




ECEC 


85J 


>l.lxlO i '' 




PHFB 


0.0525 


(a) 
(b) 


7.123x10 
1.801 xlO 22 










[59] 

EI] 
[57] 


SU(4) CTT 

MCM 

QRPA 


0.0528 
0.0071 
0.0822 


(a) 
(b) 
(a) 
(b) 
(a) 
(b) 


7.032xl0 21 
1.778xl0 22 
3.889xl0 23 
9.833xl0 23 
2.901 xlO 21 
7.336xl0 21 


1 v 
i2b Xe 


ECEC 








PHFB 


0.0487 


(a) 
(b) 


5.682xl0 24 
1.437xl0 25 


13u Ba 


(3+13+ 


[89] 
[90] 


(2.16±0.52) x 10 21t 
>4.0xl0 21 ** 


* 


PHFB 


0.0415 


(a) 
(b) 


4.797xl0 29 
1.213xl0 30 










US] 


QRPA 
QRPA 


0.0697 


(a) 
(b) 


1.700xl0 29 
4.298xl0 29 
1.370xl0 29 




(3+ EC 


OH 
M 


(2.16±0.52) x 10 2iT 
>4.0xl0 21 ** 


* 


PHFB 


0.0415 


(a) 
(b) 


4.188 xlO^ 3 
1.059xl0 24 










[59] 


SU(4) CTT 


0.0568 


(a) 
(b) 


2.235xl0 23 
5.651x10 } 










m\ 

[57] 


MCM 
QRPA 


0.0152 
0.0849 


(a) 
(b) 
(a) 
(b) 


3.121xl0 24 
7.890 xlO 24 
l.OOOxlO 23 
2.529 xlO 23 




ECEC 


ESI 
EDI 


(2.16±0.52) x 10^ iT 
>4.0xl0 21 ** 


* 

EI] 


PHFB 

SU(4) CTT 
MCM 


0.0415 
0.0568 
0.0072 


(a) 
(b) 
(a) 
(b) 
(a) 
(b) 


1.405xl0 22 
3.553xl0 22 
7.498xl0 21 
1.896xl0 22 
4.666 xlO 23 
1.180xl0 24 












QRPA 


0.0759 


(a) 
(b) 


4.199xl0 21 
1.062 xlO 22 


132 Ba 


ECEC 


[89] 

ESQ] 


(1.3±0.9) x 10 21t 
>3.0xl0 20 ** 


* 


PHFB 


0.0522 


(a) 
(b) 


5.474xl0 24 
1.384xl0 25 



* denotes present work, ** denotes (0^ + 2v) decay mode, t denotes (0^ + 2i> + QvM) decay mode 



Table 5: Effect of the variation in ( qq on (Q 2 ) , 8 2 and M 2v for 124.126,128,130^ i24,i26,i28„i30,i32 Xej and 
i30,i32 Ba nudei _ 



Nuclei ( qq 0.00 0.20 0.40 0.60 0.80 0.90 0.95 1.00 1.05 1.10 1.20 1.30 1.40 1.50 



124 Xe 


(Ql) 
Wo/ 


o 


.00 


0. 


310 


1 


.144 


36.492 


62.41 


69.16* 


71.03 


73.36 


76.17 


79.81 


85.96 


89.47*1* 


90.66 


91.87 




02 





.00 


0. 


093 





.094 


0.135 


0.190 


0.202* 


0.205 


0.210 


0.216 


0.225 


0.245 


0.259* 


0.261 


0.262 


124 Te 


(Ql) 





.00 


0. 


124 





.506 


33.151 


47.27 


52.98* 


55.03 


57.01 


58.91 


60.69 


63.98 


67.41 t 


71.15 


76.14 




Bo 

HI 





.00 


0. 


.077 


o 


.079 


0.123 


0.148 


0.157* 


0.161 


0.164 


0.167 


0.170 


0.176 


0.181* 


0.188 


0.197 




Mo 

lv± Zv 





.189 


0. 


195 





.193 


0.116 


0.075 


0.060* 


0.057 


0.052 


0.048 


0.041 


0.025 


0.023* 


0.028 


0.033 


126 Xc 


(Ql) 
Wo/ 


o 


.0 


0. 


.122 


o 


.473 


33.73 


55.53 


61.80 


64.97 


67.34 


69.26 


70.53 


77.56 


92.05 


94.74 


95.24 




02 





.00 


0. 


.088 





.092 


0.128 


0.174 


0.187 


0.194 


0.199 


0.203 


0.204 


0.218 


0.262 


0.266 


0.267 


126 Te 


(Ql) 





.00 





.058 





.541 


29.31 


41.23 


46.00 


48.04 


49.83 


51.45 


52.98 


55.32 


58.14 


61.57 


66.99 




do 
b>Z 





.00 


0. 


073 





.079 


0.114 


0.136 


0.145 


0.148 


0.151 


0.153 


0.156 


0.160 


0.164 


0.170 


0.179 




Mo 
lv± Zv 


o 


.170 





.197 





.169 


0.104 


0.077 


0.063 


0.055 


0.049 


0.047 


0.048 


0.039 


0.007 


0.002 


0.004 


128 Te 


(Ql) 
Wo/ 





.00 


0. 


.019 





.419 


2.662 


34.01 


38.08 


39.94 


41.62 


43.14 


44.61 


47.13 


49.34 


50. 81* 


52.98 









.00 


0. 


058 





.079 


0.079 


0.122 


0.129 


0.133 


0.136 


0.138 


0.141 


0.145 


0.148 


0.150* 


0.154 


128 Xe 


(Ql) 





.00 





329 





.964 


42.52 


57.17 


61.30 


62.86 


64.54 


67.90 


72.94 


90.48 


96.90 


99.83* 


100.05 




do 
HZ 


o 


.00 


0. 


091 





.092 


0.146 


0.178 


0.186 


0.189 


0.192 


0.199 


0.209 


0.257 


0.269 


0.275* 


0.275 




Mo 

lv± Zv 


o 


.142 


0. 


.138 





.139 


0.081 


0.042 


0.036 


0.033 


0.033 


0.029 


0.026 


0.003 


0.0004 


00002* 


00003 


130 Tc 


(Ql) 
Wo/ 


o 


.00 


-0.003 





.040 


1.654 


0.316 


26.80 


29.35 


32.12 


33.69 


35.03 


37.37 


39.30 


41.02 


42.65 




02 





.00 


0. 


006 





.065 


0.078 


0.079 


0.100 


0.107 


0.117 


0.120 


0.123 


0.127 


0.130 


0.133 


0.136 


130 Xc 


(Ql) 





.00 





.189 





.631 


22.26 


43.56 


48.61 


50.84 


52.17 


54.41 


55.37 


57.92 


61.13 


64.45 


77.79 






o 


.00 





091 


o 


.092 


0.103 


0.147 


0.159 


0.164 


0.166 


0.171 


0.173 


0.179 


0.186 


0.190 


2.155 




Mo 


o 


.123 





.124 





.124 


0.116 


0.061 


0.057 


0.049 


0.042 


0.038 


0.036 


0.031 


0.025 


0.025 


0.012 


130 Ba 


(Ql) 
Wo/ 


o 


.00 





.145 





.617 


56.36 


67.57 


71.39 


73.18 


79.29 


93.74 


98.74 


105.5 


108.6 


109.4 


110.5 




02 





.00 


0. 


092 





.095 


0.190 


0.212 


0.220 


0.223 


0.234 


0.272 


0.286 


0.298 


0.307 


0.310 


0.314 


130 Xc 


(Ql) 





.00 


0. 


.189 





.631 


22.26 


43.56 


47.96 


50.84 


52.17 


54.41 


55.37 


57.92 


61.13 


64.45 


77.79 




02 





.00 


0. 


091 





.092 


0.103 


0.147 


0.157 


0.164 


0.166 


0.171 


0.173 


0.179 


0.186 


0.190 


0.216 









.191 





.188 





.190 


0.028 


0.033 


0.030 


0.032 


0.041 


0.009 


0.002 


0.0002 


0.0001 


0.0002 


0.002 


132 Ba 


(Ql) 





.00 





.067 





.323 


0.750 


1.660 


33.83 


46.33 


50.29 


52.36 


54.33 


57.70 


60.30 


62.99 


64.55 




02 





.00 


0. 


085 





.094 


0.096 


0.096 


0.125 


0.163 


0.176 


0.180 


0.184 


0.191 


0.196 


0.202 


0.205 


132 Xe 


(Ql) 





.00 


0. 


.009 





.103 


0.495 


3.437 


32.64 


35.37 


37.75 


39.54 


41.12 


44.70 


46.38 


48.92 


51.33 




02 





.00 





039 





.089 


0.092 


0.092 


0.121 


0.127 


0.133 


0.137 


0.141 


0.150 


0.153 


0.158 


0.163 









.161 





.164 





.166 


0.167 


0.170 


0.142 


0.086 


0.052 


0.049 


0.046 


0.043 


0.040 


0.037 


0.034 


l50 Nd 


(Ql) 





.00 


0. 


.070 





.170 


23.55 


51.13 


64.98 


77.06 


83.75 


86.89 


87.76 


88.40 


88.93 


89.37 


89.73 




02 





.00 


0. 


066 





.075 


0.096 


0.163 


0.214 


0.257 


0.276 


0.283 


0.285 


0.288 


0.290 


0.291 


0.293 


150 Sm 


(Ql) 





.00 


0. 


104 





.649 


1.477 


33.14 


48.66 


58.45 


73.41 


83.53 


88.11 


93.23 


97.48 


99.85 


101.99 




02 





.00 


0. 


067 





.081 


0.084 


0.112 


0.155 


0.187 


0.238 


0.268 


0.281 


0.297 


0.309 


0.314 


0.319 




M 2u 





.196 


0. 


199 





.201 


0.200 


0.140 


0.113 


0.056 


0.033 


0.028 


0.027 


0.029 


0.026 


0.022 


0.015 



"denotes ( qq = 0.91, *denotes ( qq = 1.31, *denotes ( qq = 1.41. 



